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Hybrid Quarkonia with High Statistics from NRQCD 

UKQCD Collaboration - T. Manke, IT. Dmmmond, R.R. Horgan, H.P. Shanalian a * 
a DAMTP, Silver Street, Cambridge CB3 9EW, ENGLAND 

We have studied the 0(rnv 6 ) effects in NRQCD on the spectrum of heavy quarkonia and compare our results 
for different lattices (quenched and dynamical). We also report on an investigation into hybrid states within the 
framework of NRQCD. This suggests that the lowest lying hybrid is around the B* B threshold and 3 standard 
deviations above the BB. 
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1. INTRODUCTION 

NRQCD has proved to be very successful when 
applied to heavy quark spectroscopy as relativis- 
tic effects can be taken into account systemat- 
ically. In particular the low lying experimen- 
tal spectrum can be reproduced directly from 
£nrqcd- Because NRQCD allows very fast and 
accurate measurements it is also possible to ob- 
tain reliable results for the spin structure which 
reveals the quark content in the confining regime. 
In a previous paper Q we reported on a study 
of higher order corrections to the spin splittings. 
Here we will analyse these results more carefully 
on different lattices including some preliminary 
studies on unquenched configurations. It is very 
appealing to reveal also the gluon content in the 
non-perturbative regime of QCD. To this end we 
couple the magnetic field to the quark-antiquark 
pair. This hybrid state will allow quantum num- 
bers which are not permitted in the ordinary 
quark model. There are many predictions for 
the masses of such hybrids from phenomenolog- 
ical models ||J|and from lattice calculations in 
the static limit Q] . Here we go beyond the static 
approximation and employ NRQCD to treat the 
quarks dynamically. In Section 2 we present our 
Lagrangian and the results from a detailed analy- 
sis of the systematic errors in NRQCD. In Section 
3 we report on hybrid results from NRQCD. 
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2. SYSTEMATIC ERRORS 

For our investigation into the systematic errors 
of NRQCD we use the evolution equation 

G t +i = K t+1 U}(x) K t (1 - aSH) G t >i , 
G t=1 = K t=1 Ul (x) K t=0 S(x) , 
all lt 
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The source term, S(x), depends on the operator 
to be propagated and SH is given in § includ- 
ing terms up to 0(mv 6 ). The stability parameter 
n depends on the quark mass mj,. We also re- 
duced the lattice artefacts by redefining the lat- 
tice derivatives and the electromagnetic field ten- 
sor in Equation [|. All links are tadpole improved 
and we use the tree-level coefficients for all op- 
erators. For our analysis with lower accuracy in 
the velocity expansion we use the same evolution 
equation as in || . As opposed to || we employed 
gauge invariant operators throughout. We con- 
structed the meson operators with definite J PC 
from extended link variables to improve the over- 
lap with the ground state Q. 

2.1. 0(mw 4 ) vs. 0(mv 6 ) 

For T the relativistic corrections to the spin 
splittings are expected at the 10% level from 
power counting arguments. In Figure [l] we show 
our results for the hyperfine structure, 3 Si — 1 5o- 
We see a clear effect in this very accurate quan- 
tity and the change of 13(3)% compares well with 
this expectation. For charmonium the results are 
less encouraging and we see 51(2)% effects which 
are beyond the naive expectations from NRQCD. 
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Similarly, there is also a reduction of the fine 
structure, 3 Pj — P, away from the experimental 
values. This is a 34(14)% effect and it is still con- 
sistent with our expectations. But we note that 
the fine structure ratio is now consistent with ex- 
periment, Rf s = 0.56(19). This indicates that the 
inclusion of our new terms is important to repro- 
duce the internal structure of the P-wave triplet 
accurately. 3 Pj — P. 
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Figure 1. Hyperfine Structure in T. The values 
around a 2 « 0.02 fm 2 are shifted for clarity. 



2.2. (3 = 5.7 vs. f3 = 6.0 

There will be lattice artefacts from discretisa- 
tion errors which are not removed at the tree 
level. In Figure [j] and || all results are plotted 
against a 2 , determined from the IP — IS split- 
ting. We see that our new evolution equation has 
reduced the lattice spacing dependence in the spin 
structure significantly. The remaining effects can 
be interpreted as being 0(aa 2 ) |(|. A different 
tadpole prescription may reduce these discretisa- 
tion errors further Q]. Comparing our results to 
those in m we see that going from uqp ~ > u ol 
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Figure 2. Fine Structure in T. The scale is de- 
termined from the IP — IS splitting. 



changes the hyperfine splitting by the same mag- 
nitude as the change from 0(mi> 4 ) — > 0(mv 6 ) 
but in the opposite direction. We expect a simi- 
lar behaviour for the P-state triplet. 

2.3. Nf = vs. N f = 2 

We have also used preliminary data from un- 
quenched configurations at [3 = 5.2 with 2 dy- 
namical flavours 0. Our values with K sca = 0.139 
are directly comparable to the results at j3 = 
5.7, Nf = 0. For this sea quark mass we see a very 
weak dependence on the number of flavours. The 
systematic errors discussed above have certainly 
a larger effect. Therefore unqucnching effects 
should be studied for spin-independent quantities 
where spin corrections cannot obscure the analy- 
sis. However, those quantities have a bigger sta- 
tistical errors and the observed effects are small 
as well [|. 

3. HEAVY HYBRIDS 

From phenomenological models one expects the 
lowest hybrid state to come from a magnetic 
gluon J p = 1+ which couples to QQ in an octet 
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representation. The simplest lattice operators for 
such a state are 



(2) 



and the corresponding spin-singlet without er^. 
On a lattice they belong to the representations 
A~+ ,Tf+ ,T~ + ,E~+ ,T~~ . This list also in- 
cludes the exotic hybrid 1 h . As a first step in 
establishing the masses of the lowest lying hybrid 
states we employed the leading order NRQCD 
correct to 0(mv 2 ). At this level of accuracy all 
these states are degenerate in energy. The op- 
timisation of the overlap with the hybrid state 
is more elaborate and requires most of the CPU 
time. We use a fuzzing algorithm for the link 
variables from which we construct the extended 
operators with correct J PC . We calculated prop- 
agators for several combinations of operators with 
different extent in the sink and source. After av- 
eraging over 20,000 sources from 500 configura- 
tions at /3 = 6.0 we obtained a very clear signal. 
We show the correlated fit results in Figure 0. To 
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FIT RESULTS FOR HYBRID STATE 
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Figure 3. t m i n plot for a one-exponential (circles) 
and a two-exponential (crosses) fit. i max = 10. 



illustrate this result we compare the spin inde- 
pendent spectrum from NRQCD with the experi- 
mental values in Figure 0. The hybrid state is our 
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Figure 4. Spin-independent excitations in T. The 
scale is determined from the IP — 15 splitting. 



prediction and it is approximately 1.64(16) GeV 
above the T. A more careful analysis is needed 
to study the systematic errors. In particular the 
finite volume errors are of concern because such 
states will be rather broad on the lattice we used 
(16 3 x 48, a" 1 = 2.44(4) GeV). Future work will 
also include the spin-dependent terms in NRQCD 
to study the splittings within the hybrid multi- 
plets. 
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